Probing quantum and thermal noise in an interacting many-body system 
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The probabihstic character of the measurement process is one of the most puzzhng and fascinating 
aspects of quantum mechanics. In many-body systems quantum mechanical noise reveals non-local 
correlations of the underlying many-body states. Here, we provide a complete experimental anal- 
ysis of the shot-to-shot variations of interference fringe contrast for pairs of independently created 
one-dimensional Bose condensates. Analyzing different system sizes we observe the crossover from 
thermal to quantum noise, reflected in a characteristic change in the distribution functions from 
Poissonian to Gumbel-type, in excellent agreement with theoretical predictions based on the Lut- 
tinger liquid formalism. We present the first experimental observation of quasi long-range order in 
one-dimensional atomic condensates, which is a hallmark of quantum fluctuations in one-dimensional 
systems. Furthermore, our experiments constitute the first analysis of the full distribution of quan- 
tum noise in an interacting many-body system. 



I. INTRODUCTION 

The probabilistic nature of Schrodinger wave functions 
and the uncertainty principle are crucial aspects of the 
modern understanding of quantum matter. Starting with 
the famous Bohr-Einstein debates, intrinsic quantum me- 
chanical noise has been the subject of numerous discus- 
sions and controversies I*] . The analysis of quantum and 
thermal noise is not only important in the study of fun- 
damental problems in physics but also crucial for under- 
standing the ultimate limits of optical and electric detec- 
tors, sensors, and sources. 

The understanding of noise in the non-classical states 
of light facilitated the development of the theory of pho- 
todetection and led to the foundation of quantum op- 
tics [2*]. In solid state systems current fluctuations were 
used to probe the nature of electrical transport in meso- 
scopic electron systems 0, and to investigate quantum 
correlations and entanglement in electron interferometers 
[5,6]. ^ 

In atomic physics, noise correlation analysis [7[, was 
employed to study quantum states in optical lattices 
0, 9], pair correlations in fermi gases the counting 
statistics in an atom laser , and the Hanburry-Brown- 
Twiss effect for both bosons and fermions [13, [l3 • 

Interference experiments provide a different powerful 
tool which allowed, for example, the study of macro- 
scopic phase coherence [3, critical fluctuations p^ . 
thermal fluctuations in elongated condensates and 
the Berezinskii-Kosterlitz-Thouless transition in a two- 
dimensional quantum gas \vf\ . Recently it has been sug- 
gested that the full statistics of fluctuations in the con- 
trast of interference fringes can be used to probe high 
order correlation functions and reveal non-trivial phases 
of low-dimensional condensates jH, . 

In this paper we develop this idea further and provide 



a complete experimental analysis of the shot-to-shot vari- 
ations of the interference fringe contrast for pairs of in- 
dependently created one-dimcnsional Bose condensates. 
For long system sizes we find that both the average con- 
trast and its variations are dominated by thermal fluc- 
tuations. For smaller systems we demonstrate that the 
distribution functions of fringe contrast provide unam- 
biguous signatures of quantum fluctuations. Earlier ex- 
periments with one-dimensional condensates in optical 
lattices observed manifestations of strong atomic inter- 
actions in the three body recombination rate ^i]] , the to- 
tal energy (2]| , and the momentum distribution of atoms 
[2^ . However the power-law nature of the correlation 
functions, which is the hallmark of one-dimensional quasi 
condensates, has never been observed before. Our exper- 
iments provide the flrst experimental demonstration of 
quasi long-range order in one-dimensional condensates. 
This is the smoking gun signature of quantum fluctua- 
tions in systems of weakly interacting ultracold atoms. 
More importantly, our work constitutes the flrst mea- 
surements of the full distribution functions of quantum 
noise in an interacting many-body system. 



II. EXPERIMENTAL PROCEDURE 

Our experiments are performed using two indepen- 
dent one-dimensional quantum degenerate atomic Bose 

Pes created in a radio- frequency-induced micro-trap 
, on an atom chip [23 . Each sample con- 

tains typically 3000—5000 atoms in the Id regime [28. 29i, 
with both temperature T and chemical potential fulfill- 
ing ksTjfj, < hv^, where — 3.0 kHz is the trapping 
frequency of the harmonic transverse confinement. When 
the two independent Id condensates are released, they 
recombine in time-of-flight and the resulting interference 
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FIG. 1: Experimental setup and observed interference pat- 
terns, (a) Two independent Id Bose gases are created by first 
splitting a single highly elongated magnetic trap on an atom 
chip holding a thermal ensemble of atoms into a double-well 
using radio-frequency-induced potentials. In a second step the 
separate parts are evaporatively cooled to degeneracy, pro- 
ducing two individual Id condensates (schematic on the left). 
The two systems are then simultaneously released from the 
trapping potential and the resulting interference pattern is 
recorded with standard absorption imaging. The vertical ori- 
entation of the initial system is chosen so that the interference 
pattern can be imaged along its transverse direction, parallel 
to the atom chip (illustration on the right side), (b) Color 
coded images of the resulting density patterns. The observed 
interference fringes show a meandering over the length of the 
system (z-direction) , which is due to the local differences in 
relative phase between the two original Id condensates. Con- 
sequently, this waviness of the patterns contains information 
about the phase correlations in the individual condensates. 



pattern is recorded using standard absorption imaging 
(figure [TJi). 

Examples of the observed fringe patterns are shown in 
figure [Jd. While the interference patterns have high lo- 
cal contrast, the interference fringes as a whole are not 
straight lines. This meandering character of the inter- 
ference patterns shows that the relative phase between 
the two condensates is not constant but fluctuates from 
point to point. These phase variations originate from 
both quantum and thermal fluctuations in the original Id 
condensates and reflect the non-mean-field character of 
low-dimensional systems. Integrating local interference 
patterns over a finite length L leads to summing inter- 
ference patterns which are not in phase with each other 
and results in a reduction of the total fringe contrast (fig- 
ure [2]). This reduction of the interference contrast, and 
its statistical fiuctuations, contains important informa- 
tion about the phase correlations of the individual Id 
condensates and is the main quantity addressed in this 
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FIG. 2: Analysis of the observed interference patterns. For 
quantitative analysis, we integrate over central slices of vary- 
ing length L of the density profiles in the longitudinal di- 
rection as indicated by the shading in the top row, to ob- 
tain multiple transverse line density profiles. We then extract 
the interference amplitude \Aq \ by Fourier transforming these 
profiles and extracting the Fourier coefficient corresponding 
to the fringe spacing Q. To illustrate contrast reduction with 
increasing L, the fringe patterns shown in the bottom row are 
normalized. Modulated cosine fits to these profiles then yield 
contrasts C{L) which decrease with L. Note that the inter- 
ference amplitude \Aq\ as defined in equation ([T| is related 
to this contrast as \Aq\ — nidi/ x C. Consequently, as can 
be seen from equation \Aq{L)\ increases with L, whereas 
the contrast C{L) decreases with L. 

work. 



III. THEORETICAL MODEL 

Before we proceed to the statistical analysis of the ex- 
perimental data we present a quick overview of the the- 
oretical foundations of our study. The analysis is per- 
formed in two steps. In the first step we analyze the 
average amplitude of interference fringes as a function of 
the integration length L, the second step then analyzes 
the quantum and thermal noise contributions to the shot- 
to-shot fluctuations in the average contrast 

Our interference patterns show a periodic density mod- 
ulation at the interference wave vector Q = md/ht, where 
TO is the mass of the atoms, d is the in-trap separation of 
the two Id systems, and t is the expansion time. Assum- 
ing ballistic expansion, the complex amplitude of this 
density modulation after integration over a length L is 
given by [nil 

.L/2 

AQiL)= / dza\{z)a2{z). (1) 

J-L/2 

Here ai_2 are the boson annihilation operators within the 
two original one-dimensional condensates before the ex- 
pansion. The phase of Aq (L) describes the position of in- 
terference fringes and is determined by the relative phase 
between the two condensates averaged between -L/2 and 
+L/2 . 

Since we perform interference experiments with inde- 
pendent condensates, the phase of Aq is random, and 
the expectation value (Aq) is zero. This does not imply 
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the absence of fringes but shows the unpredictable ran- 
dom phase in individual interference patterns [3ll. [s^. [sst . 
Consequently, to study the contrast statistics of the 
interference patterns one has to consider the quantity 
(|Aq(L)P). This quantity is independent of the overall 
phase difference but is strongly affected by phase twisting 
within each condensate. 

In the case of ideal non-fluctuating condensates one 
expects to find perfect contrast for any size of the system. 
This imphes (|Aq(L)P) cx L^. In the opposite regime 
of short range phase correlations with finite correlation 
length the net interference pattern comes from adding 
up fringes in L/^</, uncorrelated domains. In this case 
the net interference pattern is strongly suppressed and 
appears only as a square root fluctuation, {\Aq{L)\'^) cx 

More precisely, {\Aq{L)\'^) is determined by the inte- 
gral of the two point correlation function: 

rL/2 

(|^q(L)|2) = / dzj dz2 (2) 

J-L/2 J-L/2 

{a\ (zi)ai (z2)) (4 (^2)02 {zi)). 

A special feature of one-dimensional systems of inter- 
acting bosons is the dramatic enhancement of fluctua- 
tions. Even at T = 0, true long-range order is not possi- 
ble and only quasi-condensates with a power-law decay of 
the correlation function {a'^ {z2)a{zi)) exist [33|. At finite 
temperatures one finds exponential decay of the correla- 
tion function for distances \z-i — zi\ exceeding a thermal 
correlation length ^^(r) 

To adequately describe these systems, a beyond mean- 
field theory is required. A powerful non-perturbative 
approach which describes the long distance behavior 
of the correlation functions of one-dimensional systems 
is the Luttinger liquid theory (see methods and Refs. 
[^m, [13]), on which we base our further analysis. 

Using a standard expression for the two point correla- 
tion function in Luttinger liquid theory we obtain 

{\Aq{L)\^) = L + C2 (^^^ ' '^"^ fi^,K). (3) 

Here, K w nh^J^^ is the Luttinger parameter for the 

weakly interacting Id Bose gas, with nid being the Id line 
density, g = 2hu±as the effective Id coupling constant, 
and a. being the s-wave scattering length, fh — , ^ 

is the healing length and ^^(T) = ^^'^^ is the thermal 
correlation length of the Id condensates (for the weakly 
interacting regime), while Ci and C2 are numerical con- 
stants of order unity. The function f{x, K) is given by 

.1 .1 / ^ 

f{x,K)^ / / dudv[——^\ . (4) 
Jo Jo \a;smh(-'-j — -) J 

Note that the finite number of particles can in principle 
lead to corrections to equation (jH) [s^, [s^ . We checked 



that for our parameters this shot noise is of no impor- 
tance even for the smallest L we investigate. 

The first term in equation ([3]) is generally small. It 
describes contributions from the short distance part of 
the correlation function, which is not sensitive to ther- 
mal fluctuations (provided that T < /i as in our system) 
or the non-trivial effects of quantum fluctuations in Id 
systems. The second term shows how quantum and ther- 
mal fluctuations in the Id condensates affect {\Aq{L)\'^). 

Let us first discuss the case of low temperatures and/or 
short system sizes {L/£^^{T) < 1). Here we only need 
to consider quantum fluctuations, which originate from 
interactions between atoms. Non-interacting bosons at 
zero temperature have no phase fluctuations, their in- 
terference pattern exhibit perfect contrast, leading to 
(I^q(L)P) oc L^. In the other extreme, impenetrable 
bosons (Tonks-Girardeau gas) have very strong fluctua- 
tions and their interference pattern corresponds to short 
range correlations, {\Aq{L)\'^) cx L (see discussion above, 
eq. (2) ) . For finite interaction strength we find something 
in between, resulting in the scaling {\Aq{L)\'^) (x L^~^/^. 

Finite temperature introduces thermal fluctuations, 
which create phase fluctuations with a temperature de- 
pendent correlation length ^^(T). When L > £,4,{T) ther- 
mal fluctuations dominate and the interference amplitude 
scales as {\Aq{L)\^) cx L^4,{T). 

The experimentally observed interference patterns pro- 
vide us with more information than just the average 
value {\Aq{L)\'^). As a second step of our analysis we 
consider the shot-to-shot fluctuations of individual mea- 
surements, which are characterized by the higher mo- 
ments d^Qp") and ultimately by the entire distribu- 
tion function M^(|^q(L)P). For visualizing the shot-to- 
shot fluctuations of the interference amplitude, it turns 
out to be convenient to consider the normalized variable 
a{L) = \Aq{L)\'^ / {\Aq{L)\'^) and its distribution func- 
tion W{a{L)). The importance of the higher moments 
d^Qp") is that they are directly related to the higher 
order correlation functions of the Id interacting Bose gas 

In the following, we only summarize the scaling of the 
distribution function, a formal theoretical approach is 
discussed in the methods. In the case of non-interacting 
ideal condensates, i.e. perfect interference patterns in 
each measurement, the distribution function W{a(L)) 
approaches a delta function. When interactions are weak 
but finite, we expect a narrow distribution of width 1/K 
and W(a(L)) to approach a universal Gumbel distri- 
bution |39j . In the limit of long integration lengths, 
L ^ £,<p{T), thermal fluctuations dominate. As discussed 
above, in this case the net interference pattern comes 
from adding local interference patterns from many un- 
correlated domains resulting in the distribution function 
being Poissonian. For integration lengths comparable to 
^^(T) both quantum and thermal fluctuations are impor- 
tant. In this regime we expect W{a{L)) to show a double 
peak structure, with the peak at small amplitudes coming 
from the thermal noise and the peak at finite amplitude 
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FIG. 3: Length dependence of the average contrast. The data 
points show the measured (|j4q|^) for three different temper- 
atures T. Error bars indicate s.e.m. Each data point contains 
50 individual interference measurements. The sohd hues are 
fits of equation (|3]) to the data with T as free parameter (see 
methods) and parameters Ci = 0.91 and C2 = 0.32. K and 
are determined independently from measurements of n^d 
and u±. 

from quantum noise. 



IV. ANALYSIS OF THE AVERAGE 
INTERFERENCE AMPLITUDE 

We now turn to the analysis of our experimental data, 
starting with the average interference amplitude square 
(I^Qp) ^'^d its variation with system size L. 

To extract \Aq\ from the observed interference pat- 
terns, we first obtain transverse density profiles inte- 
grated over the longitudinal direction for different length 
L, as shown in figure [2j We then fit a cosine function 
with a Gaussian envelope to the resulting fringe profiles 
to extract the relative phase and the interfering ampli- 
tude \Aq\ as functions of L (see methods). To ensure 
homogeneous Id density, we restrict our analysis to the 
central 50% of the system, where nid{z) varies at most 
by ~ 10% for the largest L considered. This modula- 
tion is neglected and we obtain a single value for ni^ by 
averaging over the atomic density in this center region. 

Figure [3] shows the experimentally observed average 
interference amplitude squares (|^qP) for three different 
temperatures, with the density nid = 50 /im^^ and the 
transverse trapping frequency v±^ = 3.0 kHz (Luttinger 
parameter K = A2 and healing length — 0.3 /im) iden- 
tical for all three data sets. The higher temperature data 
sets are obtained by waiting for different times after the 
initial preparation of the two condensates. During this 
waiting time, the system heats due to residual noise in 
the magnetic trapping fields. 

To compare measurement and theory, we fit the func- 



tion ([3]) to the experimental data (figure [3]) with the tem- 
perature T as a free parameter (see methods). We find 
the functional behavior of the measured contrasts to be 
in very good agreement with the theoretical predictions. 
This is of particular interest as the shape of these curves 
is determined by both the quantum and thermal contri- 
butions to the average contrast, as discussed above. For 
integration length longer than 20-30 /im we observe a lin- 
ear dependence of {\Aq{L)\^) on L. This corresponds to 
the L » £,,j,{T) regime where thermal fiuctuations dom- 
inate. 

For shorter segment lengths, quantum fluctuations are 
important. However the analysis presented in figure [3] is 
not sufficient to make the case for quantum fluctuations. 
The Luttinger parameter for our system is = 42, and 
it is impossible to observe the L~^/^ correction to the 
ideal case (noise-free) power law in the limited range 
of lengths available. From flgure[3]we can not prove that 
fluctuations are present at all for such short system sizes. 
We will address this in the next part of our analysis by 
demonstrating that quantum fiuctuations manifest un- 
ambiguously in the shot-to-shot fiuctuations of |Aq(L)P, 
rather than in the {\Aq{L)\'^) average value. 

^From the fits we obtain the temperatures T — 
31, 47, 65 nK for 0, 50, 100 ms waiting time, respectively. 
As expected, the greater the waiting time, the higher 
the temperature of the system. We note here that this 
method measures the temperature of collective excitation 
in the condensate. We cannot conflrm that this temper- 
ature is identical to that of the residual thermal atoms in 
the trap. Reliable detection of the thermal background 
is possible only down to T « 80 nK in our setup. 

The contrast method we present here can be used to 
measure the temperature of collective modes of Id Bose 
gases at extremely low temperatures and small atom 
numbers, suggesting the usefulness of this method for 
precise thermometry of Id condensates when conven- 
tional methods fail. 



V. ANALYSIS OF THE FULL DISTRIBUTION 
FUNCTIONS OF INTERFERENCE AMPLITUDE 

We now analyze the full information contained in the 
statistics of the interference contrast. 

In flgure IH we show histograms of the measured dis- 
tributions W{a{L)) for four different length scales L and 
two different temperatures, T = 30 nK (upper row, fig- 
ure UJa)) and T = 60 nK (lower row, figure Sljb)), ob- 
tained using the contrast average method discussed in 
the previous section. Both data sets were obtained from 
Id condensates with nid = 60 /im^^, = 3.0kHz, and 
K = AQ. 

For a detailed comparison between measurement and 
theory, we numerically calculate the distribution func- 
tions for the corresponding experiment al p arameters for 
each histogram (see methods and Ref. ^3)- We empha- 
size that once we know the temperature of the system. 
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FIG. 4: Distribution functions of the measured in- 
terference contrasts for different lengths L. (a) The 
length-dependent normalized interference contrasts a{L) — 
\Aq{L)\^ / {\Aq{L)\^) of 170 individual experimental realiza- 
tions with identical parameters [nia = 60 fim~^ , u± = 
3.0 kHz, K = 46) are displayed as histograms. The red 
curves show the corresponding calculated distributions for 
T = 30nK (C^(T) = 4.1 fim). (b) Histograms of 200 individ- 
ual measurements with the same parameters as in (a), but 
higher temperature T = 60nK {^^{T) = 2.1 urn). For both 
sets we observe very good agreement between experiment and 
theory. In particular, the predicted change of overall shape 
of the distribution functions from single peak to Poissonian 
with decreasing F = £^^{T)/L (increasing L and T) is very 
well reproduced by the experimental data. 



there are no free parameters remaining in this analysis. It 
is truly remarkable that the experimentally measured dis- 
tribution functions are in such excellent agreement with 
the theoretical prediction (figure [?]). In particular we 
clearly observe the transition from the regime dominated 
by quantum fluctuations for small L at low temperature 
to the one dominated by thermal fluctuations for large L 
and higher temperature. 

More quantitatively, the shape of the distribution func- 
tions is determined by a single dimensionless parameter 
F = K£,t/ L. For large accessed either at low temper- 
ature T or small system length L we find a single asym- 
metric peak in the distribution function W{a), resulting 
from a universal extreme-value statistics Gumbel distri- 
bution (soj which is a smoking gun signature of quantum 
fluctuations and the power-law behavior of the correla- 
tion functions. 

The Gumbel distribution has a characteristic asym- 
metric shape and typically appears when describing rare 
events such as stock market crashes or earthquakes, 
which go predominantly in one direction. This suggests 
that suppression of the contrast of interference fringes 
due to quantum fluctuations is dominated by rare but 
strong fluctuations of the phase of the bosonic flelds 
ai^2{z). Most of the time we find very small phase me- 



andering which does not affect the contrast significantly. 
Only occasionally there is a strong fluctuation leading to 
a noticeable decrease of the interference fringe contrast. 

For longer system sizes L and higher temperatures T, 
we observe that the distribution functions become Poisso- 
nian, characteristic for the dominance of thermal fluctu- 
ations and the exponentially decaying correlations. The 
system can be thought of as consisting of domains of size 
^^(r), with uncorrelated phases in each of the domains. 
In this case, adding up interference amplitudes (complex 
numbers) is similar to performing a random walk in two 
dimensions and the total amplitude (distance travelled) 
is proportional to the number of steps, resulting in the 
net interference contrast being proportional to 1/ 

Finally, for intermediate lengths and temperatures, we 
observe the formation of a double-peak structure in the 
distribution functions, characterized by a peak at zero, 
originating from thermal noise and a peak at finite am- 
plitude a originating from quantum noise jsol . [40l |. In 
this crossover regime the relative effects of quantum fluc- 
tuations are diminished but not completely suppressed, 
so that both quantum and thermal fluctuations arc of 
importance in determining the shape of the distribution 
functions. We find very good agreement between exper- 
iment and theory also in this crossover regime. 



VI. DISCUSSION AND SUMMARY 

It is interesting to note that the above analysis is based 
on the Luttinger liquid theory of interacting bosons in 
one dimension which is accurate for calculating the long- 
range part of the correlation functions but does not cap- 
ture the short distance part on the scale of the healing 
length ^/i. For system size L 3> ^/i it is this long distance 
part of the correlation functions that gives the dominant 
contribution to the integrals determining the interference 
amplitude (see equation This sensitivity to the long 
distance part of the correlation functions is the unique 
feature of interference experiments, which makes them a 
powerful tool for analyzing quantum and thermal fluctu- 
ations in low-dimensional condensates. 

Alternative approaches, such as measurements of den- 
sity fluctuations in expanding condensates [4l| are prob- 
ing correlation functions on the scale of the healing 
length. This short range part of the correlation func- 
tion is hardly sensitive to the quasi long-range nature of 
quantum fluctuations in one-dimensional systems. This 
makes it difficult to observe quantum effects by direct 
measurement of density fluctuations: they reveal the role 
of interactions but the transformation of short range 
correlations into density fluctuations masks the quantum 
correlations. 

In summary, we have studied quantum and thermal 
noise in Id systems of interacting quantum degenerate 
bosons using the full distribution function of the inter- 
ference amplitude. The shot-to-shot fluctuations in the 
contrast contain information, which can be related to 
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high-order correlation functions of the Id system. Our 
results provide the first experimental measurements of 
the full distribution function of quantum noise in an in- 
teracting many-body system. By analyzing these distri- 
bution functions we provide the first direct experimen- 
tal evidence of quasi long-range order in one-dimensional 
condensates. The remarkable agreement between our ex- 
perimental findings and theoretical predictions based on 
the Luttinger liquid model provides a confirmation of this 
theoretical approach as an effective low-energy theory of 
interacting bosons in one dimension. This demonstrates 
the power of quantum noise analysis in studying strongly 
correlated many-body systems. 

We expect our experiments to pave the way for other 
methods of characterizing many-body systems using the 
analysis of quantum noise like particle number fluctua- 
tions 43] and spin noise [3, IS, SB] ■ From the point of 
view of analyzing systems with strong interactions and 
correlations, this should allow cold-atom experiments to 
provide a complementary and different perspective to 
that provided by electron systems. 



VII. METHODS 

A. Preparing two independent Id condensates on 
an atom chip 

We start the experiment with a thermal ensemble of 

10^ Rb^"^ atoms in the \F = 2,to_f = 2 > state at a 
temperature T « 5 /iK in a single highly elongated mag- 
netic trap on an atom chip [26ll27j|. This initial sample is 
prepared using our standard procedure of laser cooling, 
magnetic trapping, and evaporative cooling [47]. The 
initial trapping configuration is then deformed along the 
transverse direction into a highly anisotropic double-well 
potential by means of radio-frequency-induced adiabatic 
potentials [23l. [25j. In particular, we employ the rf-trap 
setup presented in ref. [24] where the combination of two 
radio-frequency (rf) fields generated by wires on the atom 
chip allows the realization of a compensated symmetric 
double- well potential in the vertical plane. The final cool- 
ing of the two separated ensembles leading to the two Id 
condensates is achieved by performing forced evaporative 
cooling in the dressed state potential [4^ . We observe the 
onset of quantum degeneracy at T « 400 nK in each of 
the two potential tubes. 

The potential barrier between the two systems is con- 
trolled by the amplitude of the rf fields and the gradient 
of the static magnetic trap [1^. We realize a barrier 
height F « /cb X 4^K to ensure a complete decouphng 
of the two systems during the final cooling stage. 

After cooling and a relaxation time of 300 ms to ensure 
each system is in equilibrium (a constant rf knife is kept 
on during this time to prevent heating), each potential 
tube contains 3000 — 5000 atoms at typical temperatures 
T < 100 nK. The atoms are trapped in a strong trans- 
verse harmonic confinement of iy± 3.0 kHz at a dis- 



tance of 75 /im from the atom chip surface. Each individ- 
ual degenerate atomic ensemble is in the one-dimensional 
regime, with both temperature T and chemical potential 
fi fulfilling kBT,fi<hiy± [H,!!!]. 



B. Measuring the interference pattern 

We observe the interference pattern created by the two 
expanding, overlapping atomic clouds using standard ab- 
sorption imaging. For the vertical double-well orienta- 
tion used in the experiments, the observed interference 
fringes in the atomic density are horizontal, parallel to 
the atom chip surface. This enables us to image the in- 
terference pattern along the transverse direction of the 
system. The employed imaging system has a spatial res- 
olution of 3.4 and a noise fioor of 3 atoms per 
3x3 /im pixel. 



C. Extracting the fringe amplitude from 
interference patterns 

From a single interference image we obtain line pro- 
files for different L by integrating the two-dimensional 
absorption image over various lengths along the longitu- 
dinal direction of the system. The obtained line densi- 
ties are then Fourier transformed, and we extract Aq as 
the value of this Fourier transform at the wave vector Q 
corresponding to the observed fringe spacing. This spac- 
ing is determined from fitting the interference patterns 
with a cosine function with a Gaussian envelope plus an 
unmodulated Gaussian to account for the contrast re- 
duction. The free parameters of these fits are the rela- 
tive phase 9, the contrast, and the fringe spacing. The 
width, amplitude, and center position of the total cloud 
are determined independently from a Gaussian fit to the 
full integrated density pattern of the central area of each 
image. Note that the absolute value of the interference 
amplitude \Aq\ (as defined in eq. [1]) and the contrast C 
are related as \Aq\ = nidL x C. 



D. Average interference amplitude fits 

To compare the experimentally observed length depen- 
dence of the average interference amplitude to theory, we 
perform a least square fit of the theoretical prediction 
(equation ([3])) to the data. Since the line density nid 
is extracted directly from the absorption images and the 
transverse trapping frequency z/_l is measured precisely 
by parametric heating experiments, the only unknown 
experimental parameter is the temperature T . We in- 
clude a global additive fit parameter to account for con- 
trast reduction due to technical aspects such as finite 
imaging resolution and focal depth. 
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E. Luttinger liquid 



normalized moments of the interference fringe contrast 



A one-dimensional gas of ultra-cold bosonic atoms can 
be described by the Lieb-Liniger model of bosons inter- 
acting via a point-like repulsion [2l|, |4^, HO]- The ef- 
fective approach to the Lieb-Liniger model, capturing 
the long-distance behavior of all correlation functions is 
known as the Luttinger liquid formalism [s^ [s^ 
The essence of this approach is to represent the origi- 
nal bosonic field in terms of the two phase fields a(z) — 
{nid + 96'(z))2e*"^(^^ and keep only the terms quadratic 
in (j)[z),6{z) in the Hamiltonian. The resulting theory 
has linear spectrum of bosonic sound waves and shows 
algebraic decay of all correlation functions at zero tem- 



perature (e.g < a^a > 



Z2\^^'^^) and exponential 



decay for finite temperature (e.g. < a'(zi)a(z2) 
nid[TT / {£,TnidSmh{-n:{zi - Z2) / ir)]^^"^^ ■ Here K is the 
fundamental parameter of the theory, the so-called Lut- 
tinger parameter and — ^<j}{T)/K. The last expression 
applies down to the short distance cut-off given by the 
healing length. The value of K is uniquely determined by 
the dimensionless ratio characterizing the original micro- 
scopic model: 7 — mg/h'^nid, where g is one-dimensional 
interaction strength. In the weakly interacting regime 
studied here, K w n/y^. Recent analysis showed [51.., 5 2] 
that the Luttinger liquid formalism provides extremely 
accurate description of the correlation functions of the 
Lieb-Linger model for both long distances and distances 
just beyond the healing length. 



Calculation of the distribution functions 



(0 = (I^Qr>/(I^Qr>'"= / W{a)a"^da. (5) 

Jo 

Using Luttinger liquid theory, these moments can be ex- 
pressed [18] as the micro-canonical partition functions of 
the Coulomb gas of 2m particles 



dui...dvn 



(6) 



X exp 



Here G{x, y) is the interaction potential, which pre- 
cise form depends on the geometry of the problem 
and the temperature. At zero temperature, G{x^ y) — 
\og\x — y\, while at nonzero temperature GT{x,y) — 



log si 



sinh 



j. Real and symmetric G{x,y) can 

be decomposed as G{x,y) = J2n=T ^n'^,l{x)'i'n{y)■ 
SvLch decomposition is similar to diagonalization of a 
symmetric matrix by finding its eigenvectors and eigen- 
values. Eigenfunctions 5'„(x) and eigenvalues G„ of 
the interaction potential G{x,y) can be used to con- 
struct the height variable h{x,{tn\) = ^n^nTn(x) — 
Tn{xf/2, where r„(x) = ■^n{x)yj Gnl K, and t„ are 
fluctuating noise variables. Introducing g{{tn}) = 
J dx exp \h ( x , Hn ) II , the distribution function can be 
written as [30, kOl 



Computation of the distribution fmictions requires, in 
principle, the knowledge of all moments of the interfer- 
ence fringes amplitude. One approach to overcome this 
problem of moments was introduced in Ref. [19], where 
methods of conformal field theory and special proper- 
ties of exactly-solvable models were used to compute the 
distribution function for periodic boundary conditions at 
zero temperature. Another method, which allows to com- 
pute the distribution functions for all boundary condi- 
tions, arbitrary temperature, and in all dimensions (40| 
is based on the mapping of the problem to a generalized 
Coulomb gas model and a related problem of fluctuating 
random surfaces (for a review, see Ref. Isj). This is the 
approach which we use in our analysis. 

The full distribution function VF(q:) is defined by the 



w{a) = n 



/27r 



S[a-g{{tn})g{{-tn})]. (7) 



We compute this function using a Monte-Carlo al- 
gorithm. Random variables {tn} are chosen from 
the Gaussian ensemble, and one-dimensional integrals 
g{{tn}), gi{—tn}) are evaluated for each realization of 
{tn}- According to equation ([7]), distribution function 
W^(a) coincides with the distribution function of the 
product gi{tn})9i{—tn}) ■ In the limit of large param- 
eter F the equation ^ can be evaluated explicitly to 
show that the distribution approaches one of the extreme- 
value statistics distributions (similar to a Gumbel form) 

MM. 
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